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ABSTRACT 

We present corrections to the 



Schlegel et aL] ( [19981 SFD98) reddening maps over the Sloan Digi- 
c cap area, lb find these corrections, we employ what we dub the 



tal Sky Survey northern Galactic cap area, lb tind these corrections, we employ what we dub the 
"standard crayon" method, in which we use passively evolving galaxies as color standards by which 
to measure deviations from the reddening map. We select these passively evolving galaxies spec- 
troscopically, using limits on the Ha and [On] equivalent widths to remove all star-forming galaxies 
from the SDSS main galaxy catalog. We find that by correcting for known reddening, redshift, color- 
magnitude relation, and variation of color with environmental density, we can reduce the scatter in 
color to below 3% in the bulk of the 151,637 galaxies we select. Using these galaxies we construct 
maps of the deviation from the SFD98 reddening map at 4.5° resolution, with 1-a error of ^ 1.5 
millimagnitudes E{B — V). We find that the SFD98 maps are largely accurate with most of the map 
having deviations below 3 millimagnitudes E{B — V), though some regions do deviate from SFD98 
by as much as 50%. The maximum deviation found is 45 millimagnitudes in E{B — V), and spatial 
structure of the deviation is strongly correlated with the observed dust temperature, such that SFD98 
underpredicts reddening in regions of low dust temperature. Our maps of these deviations, as well as 
their error s, are made availa ble to the scientific community as supplemental correction to SFD98 at 
|http://www . peekandgraves2010 . coiii| 

Subject headings: astronomical databases: atlases, ISM: dust, extinction. Galaxy: local interstellar 
matter, galaxies: photometry, galaxies: elliptical and lenticular, cD, Infrared: ISM 



1. INTRODUCTION 

The Galaxy's interstellar medium (ISM) is suffused 
with particulate material called interstellar dust. This 
dust, produced by stars and in dense ISM regions, ab- 
sorbs and scatters (extincts) light and therefore contam- 
inates our detections of objects beyond the Galaxy. In 
particular, this extinction is a strong function of the fre- 
quency of the light, eliminating blue and ultraviolet light 
much more than red and infrared light. Without an ac- 
curate map of these reddening and extinction effects, ob- 
servations of distant objects have incorrectly measured 
colors and brightnesses. 

Historically, two methods were used to determine the 
extinction by Galact i c dust . Galaxy counts were used by 
Shane & Wirtanen (1967), under the assumption that 
fewer galaxies would be tound b e hind r egions with more 
Galactic dust. Knapp & Kerr (1974) used Hi column 
density, making use ot the kn own relationship b e tween 



Hi column and dust column. Burstein & Heiles ( 1978 1 
elegantly merged these two techniques to make up tor 
each of their weaknesses; the galaxy count method is 
susceptible to angular-density fluctuations caused by the 
large-scale structure of galaxies on the sky (LSS) and the 
Hi method is susceptible to variations in the dust-to-gas 
ratio. 



The release of the Schlegel et al 
dening maps represented a further 



( [19981 SFD98) red- 
and dramatic, step 



forward in our knowledge of Galactic extinction. The 
method they employed depends upon the fact that Galac- 
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tic dust radiates infrared emission and is generally opti- 
cally thin at those wavelengths. Therefore, with an in- 
frared emission map and a dust temperature map one 
can deduce the dust column, if assumptions are made 
about the radiative efficiency of the dust. For the in- 
frared maps they generated a 100 micron infrared map 
at 6.1' resolution from the combination of the Infrared 
Astronomy Satellite (IRAS) data and the DIRBE (Dif- 
fuse Infrared Background Experiment) instrument on the 
Cosmic Background Explorer ( COBE) . To determine the 
temperature they made 100 and 240 micron maps at 1.3° 
resolution generated from DIRBE. To use the notation 
of SFD98, 

E{B-V)= pIcorrX (R) , (1) 

where E {B — V) is the standard selective extinction 
measure, Icorr is the 100 micron map (corrected for resid- 
ual striation, zodical light, and the uniform extragalac- 
tic background from distant galaxies), X is a derived 
coefficient such that IcorrX represents a temperature- 
corrected infrared flux, R is the ratio of the 100 micron 
flux to the 240 micron flux, and p is a constant. 

Twelve years after the publication of SFD98, we wish 
to put these extinction measurements through a rigor- 
ous test using a principle we call the "standard crayon." 
In analogy to the standard candle principle, wherein the 
observed brightness of an object of known luminosity al- 
lows one to determine a distance, a standard crayon is 
an object of known intrinsic color whose observed color 
allows one to determine the reddening of the interven- 
ing material. We note that SFp98 has been tested in 
the past (e.g . Arce fc Goodman|l999 Chen et al. 
IStanek 1998', 'Cambresy et aL| |2(J01H Yasuda eTaT 



and Rowles fc Froebrich,2009f , but that these works typ 
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ically examine the reddening correction over small areas 
or at relatively high extinction. Our analysis is the first 
to constrain (and correct) errors in the SFD98 map over 
large areas of sky at high precision. We also acknowl- 
edge the excellent work of Schlafly et al; who have been 
conducting similar tests to our own using the blue tip of 
the stellar locus, and whose work is currently in prepa- 
ration. Our results are largely in agreement with their 
work, though each work concentrates on different aspects 
of the reddening correction and its errors. In particular, 
their data cover areas with higher extinction, and as such 
they are able to address questions of variability in the 
Galactic extinction curve. 

In this work we use passively evolving (quiescent) 
galaxies as our standard crayons, although in principle 
any object beyond the Galactic dust with a known color 
could be used to do a similar experiment. Galaxies are 
well suited to this task for several reasons. They lie out- 
side the Milky Way and thus probe the entire column 
of Galactic dust along the line of sight. They are nu- 
merous and tile the sky relatively densely. And most 
critically, their intrinsic colors can be predicted based on 
other photometric and spectroscopic galaxy properties, 
making it possible to measure deviations between their 
observed colors and their predicted intrinsic colors. 

These passively evolving galaxies, whose star formation 
has ceased, are observed to populate a tight sequence in 
the color-magnitude diagram that is quite cleanly sep- 



ies (e.g., .Strateva et al...2001j .Hogg et al. 
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20041. i<'urthermore, this "red 



sequence is tilted such that the intrinsic colors ot the 
galaxies are correlated with their absolute magnitudes. 
Thus, the magnitudes of red sequence galaxies can be 
used to determine their expected colors. Another practi- 
cal advantage of red sequence galaxies is that they have 
been surveyed in large numbers using consistent observa- 
tional methods, which minimizes the effect of systematic 
errors on relative galaxy colors. In particular, this work 



makes use o f the Sloan Digital Sky Survey (SDSS) (York 



et al. 2000 ) , which provides an enormous extant survey 
ot galaxies that have been observed in a uniform manner. 

To implement our method we must determine the 
residual variation in passive galaxy color that could be 
attributable to errors in the SFD98 reddening map. The 
first step is to move all the galaxies to a standard refer- 
ence frame using K-corrections and to remove the effect 
of extinction from dust using SFD98. After this, we de- 
termine the color variation of these galaxies with intrin- 
sic properties such as redshift, absolute magnitude, and 
galactic environment, so that we can correct for these 
variations. We also attempt to determine variations in 
galaxy color that depend on observing p aram eters, such 
as errors in the SDSS photometry (see ^3.4). After ap- 
plying these corrections, we are left with a color residual 
for each galaxy, which will be a combination of errors 
in the SFD98 reddening map, intrinsic scatter in galaxy 
color, and measurement noise. Since the intrinsic scatter 
in galaxy color dominates the effect from errors in the 
reddening map, we must average over nearby galaxies to 
determine a typical error for a region of the sky and gen- 
erate a map of errors in the SFD98 reddening correction 
map. 

In Sj2]we describe the SDSS observations and how we 



select our subsample of quiescent galaxies. In fj3]we ex- 
plain how we correct these galaxy colors to use them as 
standard crayons, as well as describe the resulting color 
distribution. In S|4]we describe our map-making method- 
ology and explore various possible sources of error in the 
maps. In j p] w e show our final maps and results. We 
conclude in ^61 

2. GALAXY SAMPLE 

There are a number of considerations that determine 
the properties of the ideal standard crayon galaxies. 

First, we require a set of galaxies that possess an in- 
trinsically tight color-magnitude relation (CMR). Sub- 
stantial scatter in the intrinsic relation has a couple of 
negative effects: it reduces our ability to accurately con- 
strain the mean color-magnitude relation, and it low- 
ers the signal-to-noise {S/N) ratio of any detected color 
residuals. This requirement has to be offset against the 
need to include enough galaxies to have reasonable statis- 
tics, as the signal we are trying to measure is weak. 

Intrinsic scatter and random noise can in theory be 
beaten down with large numbers. More problematic are 
spatially correlated color variations, since these will mas- 
querade as the variations in the foreground reddening we 
are trying to measure. There are four catagories of color 
variations that could potentially cause problems: intrin- 
sic color variations in the population, variations in the 
selection bias of the sample, variations in photometry 
errors, and variations in spectroscopic errors]^ If any 
of these four types of variations are spatially correlated, 
they will create spurious signal in the derived reddening 
maps. 

Galaxies, and red sequence galaxies in particular, are 
known to be strongly clustered in three-dimensional 
space. In addition to potentially creating spatially vary- 
ing signals that are projected onto the sky, this correla- 
tion also means that redshift-dependent variations (such 
as errors in the K-correction, or problems with detecting 
emission in the presence of strong skylines) can result in 
spatially correlated color variations. 

With these considerations in mind, this section de- 
scribes the galaxy sample selection in detail. In section 
[3] we describe further corrections to the sample in order 
to account for spurious sources of color variation that can 
contaminate our reddening maps. 

2.1. SDSS Observations 

The galaxies used in this analysis are from the SDSS 
spectroscopic Main G alaxy Survey (Strauss et al.| 2002| 
Data Release 7 (DR7, |Adelman-McCarthy eFar||2008p , 
which is an r-band magnitude-limited sample with no 
color selection. We do not use galaxies from the SDS S 



Luminous Red Galaxy sample (Eisenstein et al. 2001 1 
since these galaxies were targeted based on color cuts 
and therefore may have substantial color biases. Pho- 
tometric and spectroscopic parameters were downloaded 
from the N YU Value Added Catalog (VAGC, Blanto'^ 
et al. 20051, which uses t he updated SDSS photomet - 
ric calibration described in Padmanabhan et al. (2008). 
(We follow the colloquial terminology and refer to this 
work as "iibercal"). We note that iibercal is optimized 

^ These last are relevant because we use spe ctro scopic criteria 
to select passively-evolving galaxies (see section|2.2ll. 
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for sources with r-i color of 0.2, and our median galaxy 
r-i color is 0.4, which may introduce some very small 
systematic error. W e use the K-corrections from |Blan-| 
ton & Roweis (20071 version 4.2; earlier versions of this 
K-correction algorithm had problems in the bandpass 
projection which caused detectable systematics in our 
results. These data were supplemented with emission 
line measurements and spectral models from the MPA- 
JHU release of DR7 spectrum measurementtrl Wc also 
make use o f the proje cted 5th-nearest-ncighbor over- 



density (^s, Cooper et al. 2008) to determine the local 



environment ot each galaxy in DR7 (M. Cooper, pri- 
vate communication). Cross- matching the NYU VAGC, 
MPA/JHU spectros copic measurement catalog, and the 
( 2008 ) environment catalog produces a par- 
665,107 galaxie s that were observed as part 



Cooper et al. 

ent sample or ^ 

of the Main Galaxy Sample ( Strauss et al.|[2002 ). 



2.2. Galaxy Selection Criteria 

For this work, we wish to select a sample of passively 
evolving red sequence galaxies for our set of standard 
crayons using criteria that do not include their observed 
colors. Any selection based on galaxy color could end up 
selecting for over- or under-reddened galaxies, biasing 
the sample. Instead, we choose galaxies to be quiescent, 
non-star- forming galaxies by requiring that their SDSS 
spectra are free of detectable emission lines. 

We include only those galaxies where emission line de- 
terminations can be made: this removes galaxies without 
measured line strengths or line strength errors, low-5'/iV 
spectra (median S/N < 5 A~^), and those with redshift 
mismatches between the NYU VAGC reduction and the 
MPA-JHU reduction. This cut reduces the parent sam- 
ple to 631,949 galaxies. We examine the equivalent width 
(EW) in emission of the two strongest optical emission 
lines: Ha and the doublet [Oiil A3727, which is unresolved 
in the SDSS spectra. Figure [l] shows the Ha- [On] distri- 
bution of all 631,949 galaxies m the parent sample. There 
is a peak in the distribution near zero EW in both Ha 
and [On] ; these are the quiescent galaxies we wish to se- 
lect. There is also a narrow plume of galaxies with strong 
[On] emission but weak Ha running almost horizontally 
rightward from the zeropoint, as well as a broader plume 
extending above with substantial emission in both Ha 
and [On] . The narrow horizontal plume is dominated by 
low-ionization nuclear emission- line regions (LINERs), 
while the broad plume to the right contains star-forming 
galaxies, active galactic nuclei (AGN), and "transition 
objects," which combine both star formation and AGN 
activity ( |Yan et al |[2006] ). 

The lower panel ot I' igure [T] zooms in on the region 
around the quiescent galaxies. We assume that the peak 
of the distribution, which lies slightly off from the ori- 
gin at (0.23,0.18), represents the true zeropoint of the 
bivariate distribution, and select galaxies in an ellipse 
centered around this point (solid ellipse in Figure T]). 
The ellipse has semiaxes of 4.0 A in [On] and 0.94 A in 
Ha, chosen to enclose the majority of the quiescent ob- 
jects without substantial contamination from low-level 
star-forming or LINER-like objects, with an ellipticity 
matching the contours around (0.23,0.18). We further 
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Fig. 1. — The distribution of the 631,949 galaxies in the cross- 
matched VACG -I- MPA-JHU catalog in Ha and [On] equivalent 
width. The bottom plot is a zoom in (dashed box in top plot) on 
our selection ellipse, shown in gray in both top and bottom figures. 
The bottom plot also shows the stiffcr selection criteria (dashed 
line) that we do not use as we find them to be unnecessarily narrow. 



require that the formal errors on the Ha and [On] line 
strengths also lie within this ellipse. We experimented 
with using a more restrictive sample, represented by the 
dashed ellipse in the bottom panel of Figurejl] with semi- 
axes half the size of those of the solid ellipse, and followed 
these data through the rest of the analysis. This smaller, 
more restrictive sample produced the same results as the 
larger sample with poorer statistics. We therefore con- 
clude that the galaxy colors in the larger sample are not 
significantly more contaminated by galaxies with star- 
formation or AGN activity than the more selective sam- 
ple and proceed with the larger sample of objects defined 
by the outer ellipse in the lower panel of Figure [T] Qui- 
es cent galaxies selecte d in a similar manner were used 



Graves et al. ( 2009 ) ; figure 1 of that work shows that 



http : //wMw.mpa-garching.mpg. de/SDSS/DR7/| 



these quiescent galaxies populate the red sequence in the 
color magnitude diagram, with only a limited number of 
outliers. These criteria therefore select a relatively clean 
red sequence sample without applying an explicit color 
cut, leaving a sample of 163,101 galaxies. 
The sensitivities of the reddening corrections we mea- 
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TABLE 1 

Selection Criteria 



Quiescent galaxies 



163,101 



Good photometry in griz 159,987 
logffg < 0.7 log Fg - 1.2 
logo-r < 0.7 log Fr - 1.2 
logo-i < 0.7 log Fi - 1.15 
log cr^ < 0.6 log - 0.65 

No moving objects 155,387 
OBJC_FLAGS2 bit 3 (BAD_MOVING_FIT) not set 
OBJC_FLAGS2 bit 4 (STATIONARY) set 



No superimposed galaxies 

PetroFluXr / ModelFluXr < 1.1 

No post-starburst galaxies 



fy 



< 20% 



Final clean quiescent galaxy sample 



153,676 



151,637 



151,637 



sure in this work depend strongly on the width of the un- 
derlying color distribution of our sample galaxies. Thus, 
every effort is made to identify and exclude interlop- 
ers and color outliers, without making explicit color 
cuts that will bias the intrinsic galaxy color-distribution. 
Having identified a sample of quiescent, passively evolv- 
ing galaxies, we make several further cuts to keep the 
intrinsic color distribution small. These are listed and 
quantified in Table [l] which also lists the number of sam- 
ple galaxies remaining after each cut. In total, the cuts 
remove ~ 7% of the quiescent galaxies identified in Fig- 
ure [T] We have confirmed that the main effect of each of 
these cuts is to remove outliers from the residual color- 
color plot (see Figure |3|. 

The cuts include the following: Some imaging data 
have substantially larger-than-average photometry errors 
in the griz bands, such that they stand apart from the 
main locus of galaxies in a plot of the flux errors versus 
the measured flux. We define cuts that exclude these 
galaxies as specified in Table [ij reducing the sample size 
by ~ 2%. We further remove targets that have been 
identified as moving objects by the SDSS pipeline by re- 
quiring that the "moving fit" is good and identifies the 
object as "stationary" (another ~ 3%). Other outliers in 
residual color-color space, when examined in thumbnail 
images, have a high frequency of overlapping or super- 
imposed companion objects, which evidently affect the 
photometry. We screen out such objects by removing 
galaxies whose Petrosian (aperture-based) photometry 
differs by > 10% from the photometry based on model 
light profiles (^ 1% of the sample), as this turns out to 
be a good way of detecting overlapping or superimposed 
targets. Finally, in order to minimize the intrinsic color 
spread in our standard crayon population, we remove 
post-starburst galaxies (which may have substantially 
bluer colors) by excluding galaxies whose spectral model- 
ing includes more than a 20% flux contribution from stel- 
lar populations under 1 Gyr old (~ 1% of the sample), 
based on the "modeLcoefs" parameter in the "gaLindx" 
table of the MPA-JHU spectral catalogj^ These cuts re- 
sult in a final sample of 151,637 galaxies. 



3. STANDARD CRAYONS 

The sample described above provides a set of galax- 
ies that form a well-defined red sequence and can there- 
fore be used as standard crayons to detect reddening due 
to foreground dust in the Milky Way. There are, how- 
ever, several remaining sources of color variation that 
could produce structure in the color residual map, either 
through direct spatial correlations, or through correla- 
tions with LSS or redshift. We must identify and account 
for these variations as well as possible, then test the final 
data map for lingering spurious variations. 

These variations come in two cateogries: (1) spatially- 
dependent errors or selection biases and (2) intrinsic color 
variations within the galaxy population that correlate 
with other galaxy properties. Other sources of error are 
more difficult to correct, but we can minimize their im- 
pact in_the final residual map. We discuss these in |4.2.1| 
and |4.2.2 where we assess various sources of error in the 
residual maps. 

In this section, we consider intrinsic color variations 
due to several effects. These include the color-magnitude 
relation, environment-driven variations in galaxy colors 
at fixed mass, redshift-dependent sources of color vari- 
ation, and photometric errors. We then describe our 
method for fitting out these variations so that they do 
not bias the final residual map. 

3.1. The Color-Magnitude Correction 

The most obvious example of intrinsic color variation 
is the color-magnitude relation (OMR), in which more 
luminous galaxies have redder colors. There is also clear 
evidence that the most luminous red seque nce galaxies 
are fo und in higher-density environments (Hogg et al. 
2003| )p] Additionally there is a subtle variation of the 
CMR with redshift; galaxies at the extremes of our red- 
shift range have colors more strongly correlated with in- 
trinsic brightness. This variability in the CMR is likely 
influenced by the brightness limitations (Malmquist bias) 
of the survey. This shifting relationship may not be due 
to genuine intrinsic variations in the galaxy properties, 
but it nonetheless must be taken into account when con- 
structing our reddening maps. We note that as long as 
the variations observed in the surveyed galaxies are taken 
into account, the difference between genuine galaxy color 
evolution and the effects of observation bias will not af- 
fect our final color residuals. 

3.2. The Color-Environment Correction 

There is a well-known relationship between galaxy 
color and various measures of density of galactic environ- 
ment. There is a similar relationship within our subsam- 
ple of passive galaxies, with a slight trend of galaxies in 
dense environments being redder. To infer environ mental 
density of a galaxy, we use the 6^ parameter of |Cooper| 



5 http : //www . mpa-gEirching . mpg . de/SDSS/DR7/SDSS_indx . html 



° In fact, the distribution of typical e nvironments along the red 
sequence is complicated, showing a "saddle" shape, with ~ 
galaxies ty pically found in less dense environments than eit her 
more brTess luminous red sequence galaxies. The existence of LSS 
ther efore couples the known color variations as a function of lumi- 
nosity to the distribution of galaxies on the sky. A massive cluster, 
with a concentration of very luminous, very red galaxies would pro- 
duce the eft ect of a highly reddened sight-line if this etfect is n ot 
accounted for. 



et al. (2008) for each galaxy. The ^5 parameter is de- 
iined as the surface density of spectroscopically selected 
SDSS galaxies within 1000 km and angularly closer 
than the 5th nearest neighboring galaxy, normalized by 
redshift. We find that a simple first-order polynomial fit 
captures all the variation of color with S^. 

3.3. Method of Correction 

As will be justified in §3.5[ we fit for the variation of 
color with magnitude and environment using a median 
(minimum absolute residual) fit. Specifically we fit the 
data with 



b = ca (z) + (55/3 (z) + 7 (z) 



(2) 



where a — b represents an arbitrary K-corrected, red- 
dening corrected color, c an arbitrary K-corrected, 
reddening-corrected absolute magnitude, and a, /3 & 7 
are 8th-order polynomial functions of redshift, z, e.g. : 



(3) 



We are minimizing the sum of the absolute residuals in 
a — b, c, and 65, with the assumption of 1%, uncorre- 
lated errors i n each. This 27- parameter fit is executed 
with MPFIT ( Markwardt|2009[ an implementation of the 
Levenberg-Marquardt minimization) , first in the stan- 
dard implementation, to find a "least-squares" fit, and 
subsequently in a modified implementation to find a min- 
imum absolute residual fit. 

3.4. Photometric corrections 

The SDSS data are marvelously uniform, with excel- 
lent photometric corrections, but, at the level of our in- 
vestigation, the remaining errors in the photometry can 
be important. Once the intrinsic color relationships de- 
scribed in |3.1| and |3.2| have been removed from the 
galaxy data, our sample allows us to investigate photo- 
metric errors in SDSS that are uncorrelated to the ex- 
pected structure of fluc tuations in the SFD 98 map. 

The SDSS camera ( Gunn et al.| 1998) contains 6 
columns of CCDs, each with 5 CCDs for the 5 SDSS 
filters. Each CCD covers an area of sky 13' wide as 
the camera drifts across the sky in a series of strips, 
each dozens of degrees long. Since th e filter propertie s 



vary slightly from column to column (Doi et al. 
we compare the median of our observed galaxy 



2010[ ), 
colors 



as a function of camera column, under the assumption 
that there is no intrinsic structure in the errors in the 
SFD98 reddening map in this particular pattern on the 
sky. We indeed find that our galaxy colors do vary de- 
tectably by column from the median color, with a range 
of —0.016 magnitudes (u—z, column 2) to -1-0.027 magni- 
tudes (u— z, column 5) and a typical amplitude of 0.008 
magnitudes. We remove these spurious photometric vari- 
ations from our galaxy sample. 

3.5. Properties of the Sample 

After applying the corrections, we are left with a dis- 
tribution of residual galaxy colors. Since in principle, 
any of the 10 color pairs that can be formed from the 5 
SDSS bands can be used to determine the errors in the 
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Fig. 2. — Histograms of the error in selective extinction as de- 
termined from colors of each galaxy, shown on a log scale. One 
histogram is shown for each of the 10 SDSS colors. The Gaus- 
sian fit is overplotted, and the width of that Gaussian noted in the 
upper right of the plot. The non-Gaussianity typically begins at 
about 10% of the peak of the distribution. 



SFD98 reddening map, we are interested in the proper- 
ties of each of these residual colors. Firstly, we note that 
none of the colors have distributions that are Gaussian 
(see Figure [2]) . All distributions have wide wings that 
dominate the standard deviation, although Gaussian fits 
are relatively accurate for most colors down to about 10% 
of the peak of the histogram. We identify five separate 
sources of this non-Gaussianity, although we do not at- 
tempt to fully characterize them. Each of these sources 
can be seen in the residual g — r vs. residual r — i scatter 
plot. Figure [3j 

The first comes from interloper galaxies with residual 
populations of A stars, which form the "plume" towards 
the bottom left of the plot. These are galaxies that had a 
recent burst of star formation, but are no longer forming 
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stars. Note that this plume is dramatically reduced, but 
not removed, by the post-starburst cut implemented in 
[2.2 The second comes from galaxies with disky com- 



ponents that are viewed mostly edge on, and thus are 
significantly obscured by their own dust; this group lies 
directly along the reddening vector. The third comes 
from errors in the determination of the galaxy flux. If 
the galaxy is superimposed with another galaxy of sim- 
ilar type the color will be unaffected but brightness will 
be increased. This error will lead to an incorrect color- 
magnitude correction, which will affect the residual color. 
This is manifest as a plume toward positive g — r. This 
effect is much stronger in g — r than r — i, as the CMR has 
a much stronger dependence in g — r and thus the error in 
the correction is larger. Note that this plume i s dra mat- 
ically reduced by the photometric fit cuts in §2.2[ The 
fourth source is a population of extreme outliers, with no 
particular position in color-color space. The fifth group 
of objects that do not follow a Gaussian color distribu- 
tion, and which dominate by number, seem to simply 
follow the standard distribution of the data. It is not 
clear why the main color distribution has these "natu- 
ral" non-Gaussian wings, though they may stem from 
non-Gaussian errors in the photometry itself. 

In light of these sources of non-Gaussianity, we explore 
the use of median (least-absolute-residual) fitting. Me- 
dian fitting has the inherent disadvantage of typically 
having less overall accuracy; in a Gaussian distribution 
the error in the median is \fT^j2 ~ 1.253 times larger 
than the error in the mean. However, bootstrap resam- 
pling of the residual galaxy colors shows that the sta- 
tistical error in the median has a ~ 5% smaller scatter 
than in the mean due to the non-Gaussian distribution of 
galaxy color residuals. We could use an outlier-rejecting 
method, such as cr-clipping, to reduce the error in the 
mean. Since the color distributions seem to have intrin- 
sically non-Gaussian shapes, and it is difficult to deter- 
mine the difference between outliers and members of the 
distribution, we choose to rely on median fits throughout 
this work. 

Despite the non-Gaussian wings of the color residual 
distributions, the bulk of the galaxies do follow a Gaus- 
sian distribution, and it is these galaxies that dominate 
the uncertainty in the median. We therefore fit the color 
residual histogram with a Gaussian plus an overall ver- 
tical offset (to account for the extreme outliers) and use 
the standard deviation of the core Gaussian profile to 
calculate the uncertainty in the median. We denote the 
fitted standard deviation as tr'. 

In theory, our final reddening correction map can be 
made based on any of the measured residual colors; in 
practice it makes sense to choose the color residual with 
the narrowest distribution and therefore the lowest un- 
certainty in the median. To compare the distributions 
among colors, we use the figure of merit SE{B — V), 
the error in the selective extinction determined from the 
color of a galaxy. For a given color, b — c, corrected by 
an arbitrary magnitude band, a, we find 
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Fig. 3. — The residual color distribution of the galaxies, in g — r 
vs. r — i colors. We identify 5 contributions to non-Gaussianity 
in the distribution: internal reddening in disky, edge-on galaxies; 
superpositions of quiescent galaxies; extreme outliers; galaxies that 
contain significant A s tars; and a component that follows the main 
distribution. See 93.5|for details. 



where 5 {b — c) is the resid ual color found after apply- 
ing the corrections is ^3.3 /(a) is the ratio of extinc- 



tion in band a to E {B — V), and a (z) is the slope of 
the color-magnitude relationship, as in Equation [2] Note 
that this relationship is somewhat complicated due to 
the non-trivial projection of the reddening vector onto 
the color-magnitude relationship. Figure [2] shows the 
distribution of all ten 5E {B — V)fj_^ ^. The equiva- 
lent distributions using other bands as the reference 
magnitude (i.e., 5E [B ~ V)f^_^ g, SE {B - V),,_^^^ and 
SE {B — V)ij_^ ^, not shown) are nearly identical, though 
marginally broader. We find that all 6E [B — V) dis- 
tributions using the u band (i.e., 5E{B — V)^_^^) are 

relatively broad due to two effects: (1) the larger pho- 
tometric errors in the SDSS u band (see e.g., iibercal, 
Doi et al.|[20T0l ), and (2) the lar ger intrinsic scatter in u 
due to to variations in the stellar populations (because 
the u band is the most sensitive to small fractions of re- 
cent star formation). Colors involving the u band are 
therefore not suitable for our purposes. All else being 
equal, we would use g — r to determine reddening, as 
it has a marginally smaller fitted standard deviation in 
5E {B — V) than other colors. We note that combining 
these colors does not produce a reduced a' . 

4. EXTINCTION CORRECTION MAPS 

We wish to produce a map that covers the SDSS DR7 
contiguous spectroscopic area towards Galactic zenith 
using the residual galaxy colors we generated in ^|3j 
139,526 of the 151,637 galaxies lie in this part of the 
SDSS sky. 

4.1. Method of Map Construction 

Since we do not have a galaxy for every position in 
the sky for which we wish to compute a reddening cor- 
rection, we must apply some amount of smoothing. This 
situation is analogous to a problem that confronts single- 
dish radio telescopes, wherein a large number of telescope 
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pointings are made, recovering information only near 
that single direction. These time-ordered data (TOD) 
are generally taken in a non-homogenous pattern in the 
sky, and they must be combined to form a usable grid. 
In our case the galaxies play the role of the TOD, and 
we attempt to construct a smooth grid of reddening cor- 
rection. Typically, for any point in the grid, the values 
of nearby TODs are weighted more heavily, and distant 
TODs more weakly, with some kind of effective point 
spread function (PSF). In this work we use a Gaussian 
PSF. In contrast to the problem confronting radio tele- 
scopes, we are only collecting information from the angu- 
lar area toward the galaxy itself, rather than a large ra- 
dius around the direction in which the antenna is pointed, 
typically larger than the spacing between pointings. This 
undersampling of the sky limits the fidelity of the maps 
we construct (see [4.2). 

Since we wish to use medians, rather than means, to 
determine the reddeni ng error, we have m odified the 
single-dish code used in [Peek fc Heiies| ( pOOS] ) , SDGRID, 
to allow for weighted median tits. I'he weighted median, 
X, is the value that minimizes 



EL I 



(5) 



Wi 



Where Xi is the set of values at each point (i.e., 
6E {B — V^)f,_c a each galaxy) and Wi is the set of 
weights, in this case determined by the PSF. We have 
also modified SDGRID to determine the error in the 
weighted median at each pixel. The value of each map 
pixel is determined by a different Wi from the particular 
spatial distribution of galaxies in the vicinity. Wi varies 
due to random spatial fluctuations, edge effects, and the 
clustering of galaxies in LSS. To compute the error in 
the weighted median for each point we use a bootstrap 
method. First, we assume that the distribution of values 
contributing to each point, Xi, is relatively similar. This 
assumption allows us to replicate the measurement of x 
with a large number of sets of values Xi randomly selected 
from the data set, but using the same Wi. We are then 
able to determine the typical error in x. We note that 
our assumption of identical distributions is certainly in- 
correct (were it strictly true we would see no signal in our 
map) but that since the random galaxy-to-galaxy scat- 
ter (cr' ~ 0.03 magnitudes E{B — V)) dwarfs the map 
signal (cr' ~ 0.005 magnitudes E{B — V)), it is a reason- 
able and conservative approximation, as it only slightly 
overestimates our errors. 

Since we are constructing a map that covers a large 
area of the sky, we must take care in the choice of map- 
ping projection. Following SFD98, we use the Zenithal 
Equal Area (ZEA) projection, with a projection center 
at 6 = 90°. Since our data are all at high Galactic lati- 
tude, being projected from zenith seems a minimally dis- 
torting choice, and since the pixels in ZEA projections 
have equal area, we need not be concerned that data 
in different pixels will have systematically different error 
properties. 

In addition to choosing a PSF for our mapping rou- 
tine, we must also choose a resolution. There are two 
competing forces in our selection of resolution. At high 
resolution, we are able to capture small scale fluctuations 
in the reddening error, thus increasing the fidelity of the 



measurement. At lower resolution we have more galax- 
ies contributing to a given measurement, thus increasing 
S/N. If we assume that all spatially correlated variations 
in galaxy colors stem from the errors in Galactic red- 
dening, we can find an optimal scale with the following 
Monte Carlo test. 

We use 90% of the galaxies to build a map at some res- 
olution, O, where O is the full-width at half-maximum 
(FWHM) of the Gaussian PSF. We then take the re- 
maining 10% of the galaxies and determine the extent to 
which applying the generated map reduces the scatter: 

a' (SE [B - V),_^^^ - {SE {B - F),_,_ Je) ^ ^ 

Ce = ^ 7 ^ ^ (6) 

a' [5E{B-V\_^^^) 

Here {x)q is a map generated using the value x with a 
PSF with FWHM = 9. To produce reliable values we 
repeat the experiment 200 times and use the aggregate 
data points to compute ^e- We find that, as expected, 
increases toward small and large values of O, with a 
broad minimum between O ~ 3° and O ~ 5°. 

4.2. Sources of Error in the Map 

The beauty of our particular set of standard crayons 
is that while there could easily be complex variation in 
the colors of the galaxies, based on their metallicity, star 
formation history (see, e. g. , Figure |3]), nuclear activity, 
and internal reddening, none of these depend upon their 
position on the sky. Thus, we can ignore these variations 
and treat them simply as contributions to the random 
"cosmic" scatter, which should degrade, but not skew, 
our results. There are, however, a few concerns that can 
skew our data, where systematic errors can have an effect 
on the spatial structure on the sky. 

4.2.1. Coupling to Large Scale Structure 

In a strange role reversal, in this work the LSS man- 
ifests as an annoying systematic that gets in the way 
of our study of reddening, rather than the other way 
around. We have already dealt wit h on e effect of LSS: 
the color-density relation discussed p.2[ 

A second effect of clustering is that it couples sky posi- 
tions to specific redshifts. This correlation can be prob- 
lematic if there is variation in the color of galaxies with 
small changes in redshift, as it can contribute to erro- 
neous fluctuations in the map. We know of two possible 
causes of variations in color as a function of very small 
changes in redshift (Az ~ 0.001): K-correction effects 
and skyline effects. 

It is possible that the applied K-correction may be sys- 
tematically under- or over-correcting for the effect of the 
galaxy spectrum moving from filter to filter as a func- 
tion of redshift, contaminating our sample. Skyline ef- 
fects can also impact our galaxy colors as a function of 
small changes in redshift. The sky shines brightly with 
various molecular lines, particularly hydroxy! At red- 
shifts where either of our two selection-criteria lines ([O 
II] or Ha) align with the skylines, the errors in our deter- 
mination of those lines will dramatically increase. The 
primary effect of skyline alignment is to remove galaxies 
from our sample that we would otherwise select, because 
they do not meet the limits on EW errors in Ha. Strong 
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skylines may also allow star-forming galaxies to leak into 
the sample, if the sky is over-subtracted. To determine 
the effect of these fluctuations on our final map, we re- 
assign each galaxy the color expected from the median 
color in its Az = 0.001 redshift bin. We then rebuild the 
map with these redshift determined colors. We find that 
this map does have structure, but the variation is very 
small compared to the error in our true reddening map; 
the median absolute variations are only O.OTcr, with no 
points larger than 0.6a. We therefore can ignore the im- 
pact of fluctuations in color with small shifts in redshift. 

The third effect from LSS we consider is additional 
unmodeled noise in the map. To illustrate this effect, we 
propose the following situation. Galaxies that contribute 
weight to map pixel A are more clustered on the sky than 
those that contribute to map pixel B. Imagine the true 
sky reddening has significant fluctuation from the SFD98 
reddening map below the resolution of our map. This 
sub-resolution fluctuation will add error to both pixels, 
but since the errors are significantly more correlated from 
galaxy to galaxy in pixel A, the error will be higher in 
pixel A. To tes t the effects of this, we make a map, using 
the method of |4.1[ of the angular distance to the nearest 
neighboring galaxy 4>min- We then examine the reduced 
of our reddening map as a function of this "clustering 
map", {(j>min)0- We find that, even for small values of 6 
where we expect this effect to be strongest, there is no 
evidence for increased noise with increased clustering. 

4.2.2. Large-scale striation from systematic photometric 

errors 

Ubercal made great strides in reducing and character- 
izing the photometric errors in the SDSS data. While 
iibercal reduced the errors to less than 1% in g,r,i, and 
z bands, it also uncovered very small systematic varia- 
tions that depended upon the stripes in the data that it 
could not correct. SDSS photometry data are taken in 
long tracks, and iibercal found that while it did signif- 
icantly reduce the variation in photometric values from 
track to track, detectable striations remained (see §5.4 
and figure 13 of iibercal), especially in z band. To ex- 
amine the effect of these striations in our data, we make 
maps at a variety of resolutions in SDSS A and 77 co- 
ordinates, where A is along the scan direction and 77 is 
across the scan direction (see Figure 4|. We then take 
the largest fully sampled square of data, and examine 
the Fourier power in 2D, using the method described in 
SFD98, §3.1.1 to avoid edge artifacts (Figure III). Con- 
sistent with the findings of iibercal we find thafthe only 
detectable striations occur along the A axis and that col- 
ors involving z are significantly more striated that those 
involving only g, r, and i. In particular we find that all 
maps have detectable striation at = 1°, no maps have 
detectable striation at Q = 10°, and only the g — r maps 
lack stria tion within the optimum range 3° < < 5° 
found in §4.1| We therefore adopt 8 = 4.5° as our map 
resolution and g — r as the color by which we measure 
E{B-V). 

5. STRUCTURE OF THE EXTINCTION CORRECTION 

MAPS 

Our final dust correction maps (and associated errors) 
are shown in Figure [Hj and are made available to the 
public at http : //www . peekaiidgraves2010 . com. We de- 



tect variations from the SFD98 reddening map at greater 
than 3f7 over large areas of the map, with detection 
greater than lOfi in regions toward the Galactic anti- 
center at lower latitudes. The map has Icr sensitivity of 
1.5 millimagnitudes in E (B — V), with maximum devi- 
ation of 45 millimagnitudes toward a — 128°, 5 — -f62°. 
We note that this position is coincident with both the 
ultra-faint dwarf Ursa Major II and a bright knot in the 
High- Velocity Cloud (HVC) complex A. While it is pos- 
sible that this galaxy or HVC complex is contributing 
to error in the Galactic reddening, we find no associa- 
tion between our residual reddening map and either other 
ultra-faint dwarfs or other HVC complexes. We therefore 
expect that this association is simply a coincidence. 

While deviations from SFD98 are small in absolute 
scale, they can be fractionally large, exceeding 50% in 
excess or decrement in some areas. On the whole, the 
map also shows that towards lower latitudes (and higher 
extinction) the SFD98 map underpredicts extinction rel- 
ative to higher latitudes. This trend toward underpre- 
diction of reddening in regions of increased dust is not 
in disagreement with the claims that SFD98 overpre 
diets the dust in high extincti on regions (e.g. 
Goodman 1999] |Stane k| [T998l ICambresy et al 
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we do not have data in regions with extinction above 
E {B — V) = 0.15, and these works typically examine 
areas with E {B — V) > 0.5. We do not find obvious cor- 
relations of the map with the IRAS 1 00 um map . Galac- 
tic Hi column density maps ( Kalberla et al.|2005 1 , or the 
LSS in the SDSS. 

It is also worth noting that the application of our cor- 
rection to the SFD98 maps results in only 40 square ar- 
cminutes (^7 SFD98 map pixels) of area with a negative 
reddening correction, the largest of which is only —0.3 
millimagnitudes E(B — V), far smaller than our error 
bars. For these very few areas we report a final redden- 
ing correction of zero. This result can be taken as at 
least a partial validation of the SFD98 zero-point; were 
the zero-point set incorrectly too low we would expect to 
find large areas with negative reddening. 

5.1. Correlation With Dust Temperature Map 

For comparison to our reddening correction maps, we 
generate a map of the dust temperature. This map is 
generated by determining the dust temperature at each 
galaxy position using the SFD98 temperature map, as 
derived from the DIRBE instrument data. We then take 
these values and apply the identical map construction al- 
gorithm described in |4.1[ This map is shown in Figure 
W\ We note that this map, with red representing cooler 
dust, is very similar map shown in FigurelO] with features 
reproduced down to the resolution scale. This correlation 
implies that the SFD98 map underpredicts reddening in 
regions of low dust temperature. In the context of Equa- 
tion [l] X (R) is too small for regions of relatively low R. 
Wc note that this similarity does not persist through- 
out the entire map, diverging most prominently toward 
Galactic center. Perhaps this divergence indicates that 
the X (R) function used is more appropriate for that area 
of sky than toward other areas. It is also possible that 
the superposition of grains of different temperatures (see 
figure 2 of SFD98) or varying grain morphology or com- 
position could be the culprit. The correlation between 
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our residual map and the DIRBE temperature map is 
also strong evidence for the claim that we are indeed 
measuring errors in the reddening map, rather than an 
unrelated systematic. 

6. CONCLUSIONS 

In this work we used quiescent galaxies as "standard 
crayons" , and found errors in the SFD98 reddening map 
at high Galactic latitude. We found that a pure spectro- 
scopic selection of quiescent galaxies with no Ha or [On] 
emission produced a very homogeneous population of 
galaxies. Once these galaxies were corrected for redshift- 
dependent color-magnitude relations and color-density 
relations, we found they had a scatter of less than 3% 
in the bulk of galaxy colors. We determined that since 
there were significant departures from an intrinsic Gaus- 
sian color scatter, a median fitting method was desir- 
able over a standard fit. We found that systematic er- 
rors that stem from large-scale structure in the universe 
did not impact our data, and that our data were rela- 
tively robust to striation from systematic errors in the 
SDSS photometry at a resolution of 4.5° and using g — r 
colors. Our final correction maps showed that SFD98 
maps have errors ranging from 45 to —10 millimagnitudcs 
in E{B — V), although the bulk of the map has errors 
less than 3 millimagnitudes. These maps are provided 
to the public at http://www.peekandgraves2010.com, 
Wc find the largest errors in the SFD98 map arc towards 
lower latitudes. We find some noticeable correlation with 
the temperature map SFD98 derived from DIRBE data, 
particularly toward the outer galaxy, suggesting that the 



function relating temperature to emission correction fac- 
tor employed in SFD98 was in error for dust grains to- 
ward the outer Galaxy at high latitudes. 



Funding for the SDSS and SDSS-II has been pro- 
vided by the Alfred P. Sloan Foundation, National Sci- 
ence Foundation, the U.S. Department of Energy, the 
National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Council for England. The 
SDSS Web Site is |http : / /w ww . sdss . org/ , 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Par- 
ticipating Institutions are the American Museum of Nat- 
ural History, Astrophysical Institute Potsdam, Univer- 
sity of Basel, University of Cambridge, Case Western 
Reserve Drexel University, Fermilab, the Institute for 
Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and 
Cosmology, the Academy of Sciences (LAMOST), Los 
Alamos National Laboratory, the Max-Planck-Institute 
for Astronomy (MPIA) , the Max-Planck-Institute for As- 
trophysics (MPA), New Mexico State University, Ohio 
State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 

JEGP wishes to thank Doug Finkbeiner and Eddie 
Schlafiy for very hepful conversations regarding both 
SDSS and SFD98 and for their hospitality. 



REFERENCES 



Adelman-McCarthy, J. K., Agiieros, M. A., AUam, S. S., AUende 
Prieto, C, Anderson, K. S. J., Anderson, S. P., Annis, J., 
Bahcall, N. A., Bailer-Jones, C. A. L., Baldry, I. K., Barentine, 
J. C, Bassett, B. A., Becker, A. C, Beers, T. C, Bell, E. P., 
Berlind, A. A., Bernard!, M., Blanton, M. R., Bochanski, J. J., 
Boroski, W. N., Brinchmann, J., Brinkmann, J., Brunner, 
R. J., Budavari, T., Carliles, S., Carr, M. A., Castander, P. J., 
Cinabro, D., Cool, R. J., Covey, K. R., Csabai, I., Cunha, 

C. E., Davenport, J. R. A., Dilday, B., Doi, M., Eisenstein, 

D. J., Evans, M. L., Pan, X., Pinkbeiner, D. P., Priedman, 

S. D., Prieman, J. A., Pukugita, M., Gansicke, B. T., Gates, E., 
Gillespie, B., Glazebrook, K., Gray, J., Grebel, E. K., Gunn, 
J. E., Gurbani, V. K., Hall, P. B., Harding, P., Harvanek, M., 
Hawley, S. L., Hayes, J., Heckman, T. M., Hendry, J. S., 
Hindsley, R. B., Hirata, C. M., Hogan, C. J., Hogg, D. W., 
Hyde, J. B., Ichikawa, S.-i., Ivezic, Z., Jester, S., Johnson, 
J. A., Jorgensen, A. M., Juric, M., Kent, S. M., Kessler, R., 
Kleinman, S. J., Knapp, G. R., Kron, R. G., Krzesinski, J., 
Kuropatkin, N., Lamb, D. Q., Lampeitl, H., Lebedeva, S., Lee, 
Y. S., Leger, R. P., Lepine, S., Lima, M., Lin, H., Long, D. C, 
Loomis, C. P., Loveday, J., Lupton, R. H., Malanushenko, O., 
Malanushenko, V., Mandelbaum, R., Margon, B., Marriner, 
J. P., Martmez-Delgado, D., Matsubara, T., McGehee, P. M., 
McKay, T. A., Meiksin, A., Morrison, H. L., Munn, J. A., 
Nakajima, R., Neilsen, Jr., E. H., Newberg, H. J., Nichol, 
R. C, Nicinski, T., Nieto-Santisteban, M., Nitta, A., Okamura, 
S., Owen, R., Oyaizu, H., Padmanabhan, N., Pan, K., Park, C, 
Peoples, J. J., Pier, J. R., Pope, A. C, Purger, N., Raddick, 
M. J., Re Piorentin, P., Richards, G. T., Richmond, M. W., 
Riess, A. G., Rix, H.-W., Rockosi, C. M., Sako, M., Schlegel, 
D. J., Schneider, D. P., Schreiber, M. R., Schwope, A. D., 
Seljak, U., Sesar, B., Sheldon, E., Shimasaku, K., Sivarani, T., 
Smith, J. A., Snedden, S. A., Steinmetz, M., Strauss, M. A., 
SubbaRao, M., Suto, Y., Szalay, A. S., Szapudi, L, Szkody, P., 
Tegmark, M., Thakar, A. R., Tremonti, C. A., Tucker, D. L., 
Uomoto, A., Vanden Berk, D. E., Vandenberg, J., Vidrih, S., 
Vogeley, M. S., Voges, W., Vogt, N. P., Wadadekar, Y., 
Weinberg, D. H., West, A. A., White, S. D. M., Wilhite, B. C, 
Yanny, B., Yocum, D. R., York, D. G., Zehavi, L, & Zucker, 
D. B. 2008, ApJS, 175, 297 



Arce, H. G., & Goodman, A. A. 1999, The Astrophysical Journal, 
512, L135 

Baldry, L K., Glazebrook, K., Brinkmann, J., Ivezic, Z., Lupton, 
R. H., Nichol, R. C, & Szalay, A. S. 2004, ApJ, 600, 681 

Balogh, M. L., Baldry, I. K., Nichol, R., Miller, C, Bower, R., & 
Glazebrook, K. 2004, ApJ, 615, LlOl 

Blanton, M. R., & Roweis, S. 2007, The Astronomical Journal, 
133, 734 

Blanton, M. R., Schlegel, D. J., Strauss, M. A., Brinkmann, J., 
Pinkbeiner, D., Pukugita, M., Gunn, J. E., Hogg, D. W., Ivezic, 
Z., Knapp, G. R., Lupton, R. H., Munn, J. A., Schneider, 
D. P., Tegmark, M., & Zehavi, I. 2005, AJ, 129, 2562 

Burstein, D., & Heiles, C. 1978, Astrophysical Journal, 225, 40, 
a&AA ID. AAA022. 131.051 

Cambresy, L., Boulanger, P., Lagache, G., & Stepnik, B. 2001, 
Astronomy and Astrophysics, 375, 999 

Chen, B., Pigueras, P., Torra, J., Jordi, C, Luri, X., & 

Galadi'-Enri'quez, D. 1999, Astronomy and Astrophysics, 352, 
459 

Cooper, M., Newman, J., Weiner, B., Yan, R., Willmer, C, 

Bundy, K., Coil, A., Conselice, C, Davis, M., Paber, S., Gerke, 
B., Guhathakurta, P., Koo, D., & Noeske, K. 2008, Monthly 
Notices of the Royal Astronomical Society, 383, 1058 

Doi, M., Tanaka, M., Pukugita, M., Gunn, J. P., Yasuda, N., 
Ivezic, Z., Brinkmann, J., dc Haars, E., Kleinman, S. J., 
Krzesinski, J., & Leger, R. F. 2010, arXiv, astro-ph.IM 

Eisenstein, D. J., Annis, J., Gunn, J. E., Szalay, A. S., Connolly, 
A. J., Nichol, R. C, Bahcall, N. A., Bernardi, M., Buries, S., 
Castander, P. J., Pukugita, M., Hogg, D. W., Ivezic, Z., Knapp, 
G. R., Lupton, R. H., Narayanan, V., Postman, M., Reichart, 
D. P., Richmond, M., Schneider, D. P., Schlegel, D. J., Strauss, 
M. A., SubbaRao, M., Tucker, D. L., Vanden Berk, D., Vogeley, 
M. S., Weinberg, D. H., & Yanny, B. 2001, AJ, 122, 2267 

Graves, G. J., Faber, S. M., & Schiavon, R. P. 2009, ApJ, 693, 486 



10 



g-r g-i g-z r-z 




■30 -20 -10 10 20 30 -20 -10 10 aC 3D -20 -10 10 20 30 -20 -10 10 20 30 

A [deg] 

Fig. 4. — Map of the redding in SDSS coordinates rj and A over the largest contiguous square area of the SDSS Northern Galactic cap. 
We show reddening as predicted from the 4 color pairs that have tr' < 0.03 and for 3 different values of ©. Note the striation in redder 
color pairs and for lower values of 0. 
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Fig. 6. — Our final extinction correction map. We show both our extinction correction and our significance contours. Galactic longitude 
and latitude are overlaid, with Galactic zenith at the middle of the map and Galactic center off the right side of the map. Note that the 
sign of this map is such that it portrays the residual errors in the SFD98 map, or that which must be subtracted to improve the SFD98 
map. 




Fig. 7. — Temperature map of SDSS DR7 contiguous area, generated by finding the temperature value in the DIRBE map toward each of 
our galaxies and then applying the same method as used for our reddening correction map to facilitate comparison. Red represents cooler 
areas. Note the similarity to Figure ro] in regions away from Galactic center (left side of the map). 
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